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I. INTRODUCTION 

Starting from the pioneering paper of H. Blum [l3[ . the 
medial axis as a descriptor and classifier of shapes and 
figures has been established as the best defined and stud- 
ied mathematical concept in reference to thinning and 
skeletonization of contours and shapes |E S S S Uli EH 
[T^ . [Til l . From the various mathematical tools (Maximal 
disks, cut loci, Voronoi sets 0, 0, H, H, E3, E3, [3 E3, , 
the maximal disk method seems to be the most well 
studied and applied, both in mathematical definition and 
properties [EH, El EI an< ^ m applications 
Its definition is best presented in the following form: 

Definition 1.1. Let D be a subset of ffl 1 . A closed ball 
B (or disk in 2-D) is said to be maximal in D if it is 
contained in D and if B C B' , where B' is another closed 
ball also contained in D, then B = B' 
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=(x,y) = max{\x! - x 2 \, \yi - y 2 \} 



(2) 



Called here the maximum coordinate metric. And the 
addition metric: 



rfi(x,y) = {\xi - x 2 \ + \yi - y 2 \} 



(3) 



The above definition cannot be applied to one of the 
above given metrics. The results would not lead to a 
proper skeleton. In most cases it would lead to discon- 
nected medial axes, contrary to the definition of a skele- 
ton. On the other hand, the Euclidean matrix has seri- 
ous problems when applied in discrete bitmap images [9] . 
Furthermore, due to its use of a square root function, it 
is more calculation intensive than the other metrics. 

In a completely different field and with entirely differ- 
ent motivations, some similar concepts were developed 
by P. C. Stavrinos in the early 80's 0- Q The author 
was attempting to develop better tools for the classifi- 
cation and extraction of features of various geometric 
constructions, such as classes of two dimensional man- 
ifolds immersed in a three dimensional Euclidean space. 
The tools would be developed for applications in vari- 
ous branches of mathematics and physics, for example 
in knot theory, convexity, flows, the study of differen- 
tial equations and the propagation of their solutions and 
correspon ding singularities, as described by the Huygcns 
principles [TEL fl6j . 

The developed concepts described global characteris- 
tics of surfaces connecting them to local geometric and 
topological features. The concepts presented in 0, EJ, 
mainly what the author called the first curvature o 
cally convex type, were extended and generalized in 
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and some new concepts, presented here in Section 2, were 
introduced there for the first time. At that time, the au- 
thors of and p| were largely unaware of the potential 
of their work in morphological applications. Their inter- 
est was mainly in developing new mathematical concepts 
of a not quite local but not fully global nature, as tools 
in the study and classification of geometric objects and 
properties. Some such results on non-analytic curves em- 
bedded in 5ft ra , convexity and knots' theory will be pub- 
lished elsewhere. 

When the research was turned towards morphology, 
the concept of the contact disk and the Divider set, as an 
alternative skeleton method, were presented in The 
concepts were presented in a discrete lattice environment 
and the maximum coordinate metric was used. 

In this work, a full definition of the concepts of the 
contact curvature, the contact circle and the Divider set 
will be presented. Emphasis will be put on the Euclidean 
metric method, on the Euclidean plane and three dimen- 
sional space, although, in the discussion sector, a compar- 
ison with the discrete version will be attempted, with ref- 
erence to specific applications, such as OCR and robotic 
navigation. In Section^ a brief review of previous ideas 
and definitions will be given, based mainly on In 
Section UTTI the basic definitions will be presented, in the 
form of a set of non linear algebraic equations and in- 
equalities. Some fundamental results will be presented. 
Discussion and conclusions will be the subject of Section 

hyi 

II. THE CURVATURE OF LOCALLY CONVEX 
TYPE 

In previous work 0, Q , the following definitions were 
given: 

Definition II. 1. A surface S in the Euclidean space 5ft 3 
will be of locally convex type iff, Vp G S, given < ro 
, and r G (0,ro) the intersection B 3 (p,r)f s \S, where 
B 3 (p,r) a closed ball with center p and radius r, will be 
the closure of a simply connected coordinate neighborhood 
of p on S. 

Definition II. 2. Vp G 5, the curvature of locally convex 
type, ffi 1 ' will be defined as the inverse of the infimum 
ri.c.t. of all radii of closed balls B 3 (p,r) for which the 
intersection B 3 (p, r) P| 5 will be disconnected^] . 

In the above, 5 was a connected, piece by piece smooth 
surface, which means that it was the union of a finite set 
once continuously differentiable surfaces 0, ■ 

In , an extension of the curvature of locally convex 
type was given. The definition now included any point 
p G 5ft 3 : 

Definition II. 3 (a). Vp G 5ft 3 , the curvature of locally 
convex type, Ki, c .t, will be defined as the inverse of the 
infimum ri. c .t. °f oil radii r of closed balls B^ 3 \p,r) for 
which the intersection B^ 3 > (p, r) p| 5 will be disconnected. 



ri.c.t = infim{r 1 < r — > B^(p, r) (^j 5 disconnected} . 

Another generalization of Def. 2.1 is now possible: 

Definition II. 4 (b). A surface S in 5ft 3 will be of locally 
convex type if there is no sequence of points pi G 5ft 3 , 
i £ K, the set of natural numbers, converging to a limit 
point pgSR 3 , where < K (p) < oo, such that: 

lim K{pi) < oo (4) 

This definition should be compared with assumptions 
about the suitability of analytic curves, as stated in 0/. 
Whether this definition is equivalent to analyticity, isfor 
now an open question for the authors of this work. Jj/ 

The next result has to do with the set of all points in 
5ft 3 for which < Ki. c £. This set has been named 11(5). 

Theorem II.5. The set 11(5) of a curve S in 5ft 2 or 5ft 3 , 
or a surface S in 5ft 3 is open by the Euclidean topology^. 

In Q, a series of results follows, concerning necessary 
and sufficient conditions for a point p G 5ft 2 or 5ft 3 to be- 
long to 11(5), 5* being a curve or a surface. The most 
important and utile results will be presented here, with- 
out proof, as follows (For detailed proofs, see p|): 

Lemma II. 6. If and only if a curve S is (a) a circle S 1 , 
or (b) a connected one dimensional subset of a straight 
line, the set 11(5) will be empty. 

Lemma II. 7. If S is a curve in 5ft 2 and is homeomorphic 
to (0,1), (0,1] or [0,1], then a necessary and sufficient 
condition for a point p G 5ft 2 to belong to 11(5) is that a 
compact, connected, proper subset Q of S exists, which is 
a circular arc with p as its center and disconnects 5 by 
not containing any endpoints. Therefore the relation: 

d(p,q) = constant holds V q G Q{d(p, q)} 

Furthermore, there exists an open neighborhood N(Q) 
of Q such that \fq G Q and q' G (N(Q) — Q), it is true 
that d(p, q 1 ) < d(p, q). 

The compact subset Q may be a single point, which is 
the case if 5 is analytic and contains no constant curva- 
ture arcs, as postulated in for example. 

Lemma II. 8. If S is a curve in 5ft 2 and is homeomorphic 
to S^' , then a necessary and sufficient condition for any 
point p G 5ft 2 to belong to 11(5) is either the existence 
of two points q\, qi, and corresponding open neighbor- 
hoods Ni(qi) and A^fft) such that: Vg,' G (-/Vi(Qi) — Qi), 
the relation: d{p 1 q' i ) < d(p,q) will hold, i = 1,2, or 
the existence of two compact, connected, proper subsets 
Qi, with corresponding open neighborhoods Ni(Qi) will 
exist, which will be circular arcs having p as a center 
and the relation of Lemma \II.7\ holds: Vg; G Qi and 
Vg- G (N(Qi) - Qi), it is true that d(p, g-) < d(p, q^, 
i = l,2. 
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The conclusion of Q (Theorem 3.11, p. 279), was 
that, for the calculation of the curvature of locally con- 
vex type, JsTi.c.t of a curve S, at any point p € 5ft 2 , all 
perpendiculars from p to S should be found and those 
for which the following property holds, should be se- 
lected: d(p,ql) > d[p,qi), where qi, q[ as in the previous 
lemma, t = l,2,...,n, if n such perpendiculars exist. 
Then the distances should be put in increasing order, 
d(p,qi) < d(p,q2) < .... Then, the curvature of locally 
convex type would be the inverse of the second distance 
in the ascending order series: 

K Lc .t{p) = , 1 . 

d{p 1 q2) 

A simple corollary of the above results is the close cor- 
respondence of Tl{S) and the curvature of locally convex 
type to the evolute of a curve. For example, the ellipse 
has as 11(5') the area enclosed by its evolute ITTT1 A simi- 
lar result is true of the parabola IIIII Specifically, the set 
II(S'), if 5* is a parabola, is the area bounded by the con- 
vex side of the evolute. As it will be demonstrated in the 
next Section IlIII the symmetries of a curve, the symme- 
tries of its evolute and some other geometric characteris- 
tics will decide the topological properties of the Divider. 
A very characteristic case is the hypotrochoid and its 
evolute HTT1 In this case, the symmetry of the curve and 
its evolute and most of all, the multiplicity of branches of 
the evolute which turn their convex side to a certain area 
will decide the multiplicity of branches of the Divider 
radiating from the center of the hypotrochoid 4] IIIII 

The above results were applied to the derivation of the 
concepts presented in the next Section IIIII These are 
the contact directions, the contact circles, the contact 
curvature and the Divider set. 



III. THE CONTACT DISKS AND THE 
DIVIDER SET OF A CURVE OR SURFACE 

A series if definitions is in order. They are based on 
the previous concepts and are a natural continuation of 
the previous theory. 

Definition III.l. Let a curve S in 5ft 2 be once contin- 
uously differentiable and have a perpendicular at each 
point. Let p 6 S and let B(k,r) be a closed disk 
of center k and radius r = d(k,p). Furthermore, let 
Sf]B(k,r) = {p}. Then the direction \p,k\ will be called 
a contact direction of S in p. 

A similar concept can be defined for a curve or a surface 
S in 5ft 3 . Then the two dimensional disk is replaced by a 
three dimensional ball. 

The fact is that in every case the point p is the only 
common point of S and B(k,r). lip is not an end point 
and given the assumptions of a continuous first deriva- 
tive and an unambiguously defined perpendicular, there 
are two contact directions at each p in a curve in 5ft 2 or a 



surface S in 5ft 3 . They are the two directions of the per- 
pendicular. If, on the other hand, S is a curve in 5ft 3 , then 
there is a plane normal to the curve and every direction 
on this plane is a contact direction. 

But there are extreme cases, such as end points or sin- 
gularities. In those cases, there will be a whole cone of 
contact directions at p. If p is the end point of a plane 
curve, there is a semicircle of directions with its diameter 
on the perpendicular and at the side of it pointing away 
from the curve. If the curve or surface lies on 5ft 3 , then 
there is a hemisphere of contact directions. In general, 
the space of contact directions depends on the morphol- 
ogy of the curve or surface S in the neighborhood of 
p. The same goes for singularities, where there is not a 
uniquely defined tangent and perpendicular at p. If at 
a point p on 5* a convex angle exists on one side and, 
naturally, a concave angle exists on the other side of S, 
then the contact disk on the convex side will have infi- 
nite curvature and the contact directions on the other 
side will cover an area, disk section in 5ft 2 or ball section 
in 5ft 3 , whose shape will depend on the morphology of the 
singularity (See the definitions of sharp and dull corner 
points in 

In each and every one of the above cases, the set of all 
disks or balls, having a single contact point p with S, is 
well defined. If the disk is enlarged, by moving its cen- 
ter k away from p and thereby enlarging its' radius \p,q\, 
one of two things will happen. Either, at some point, 
the intersection Sf]B(k,r) will contain other points of 
S, besides p, or the intersection will contain only p all 
the way until k reaches infinity. In the former case, there 
must be a supremum of the radii of such disks having the 
property Sf]B(k,r) = p. In the second case the supre- 
mum is infinite. Then we have the following definition. 

Definition III. 2. The inverse of the radius of the supre- 
mum of all disks which have the property: 

Sf)B(k,r) = {p} 

is called the contact curvature of S at p. The locus 
of all centers of the supremum disks, Vp £ S, is called 
the Divider set of S. 

It is obvious that the set of all disks with the property 
Sf]B(k,r) = {p} is linearly ordered by inclusion. This 
means that, as the radius of each disk increases, it con- 
tains all disks with smaller radii. It is equally obvious 
that, if S is a closed curve, in the sense of being homc- 
omorphic to containing an area, and the contact 

direction is taken towards the inside of the area, then all 
disks with the above property are contained in the area 
surrounded by S. Then, assuming that the supremum 
is not a maximum, the above definition is identical with 
that of a maximal disk [(J 0, H, 0, El El ■ Also obviously, 
since the contact direction is defined for both sides of a 
closed contour, the Divider is equally well defined for an 
arc homcomorphic to a one dimensional connected sub- 
set of a straight line or a one dimensional circle it 
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may cover cases where the maximal disk definition can- 
not be applied easily, or not at all. Furthermore, since 
the above definitions may easily be extended to curves 
with multiple cut points, II I II or even families of discon- 
nected curves or surfaces, the Divider concept can cover 
all relative medial axis definitions, such as the Voronoi 
set, etc. This advantage is expected to be decisive in 
the case of the morphological study of a whole page of 
handwritten text, as a first step for handwriting OCR. 

By far the most important advantage of the Divider 
definition is the fact that it can be easily and naturally 
extended in a discrete lattice environment, by the use 
of another metric than the Euclidean one jfj • The rules 
for the construction of the contact disks and the Divider 
in [j| are a first tentative attempt for the development 
of algorithms suitable for the discrete lattice case. The 
advanta ges a nd potential applications of the discrete lat- 
tice [Tl|, [l2j, will be included in these authors' future 
work. As already mentioned in the applications con- 
sidered so far are OCR on handwritten text in a dis- 
crete lattice environment and robotic navigation and map 
making in complex enclosed spaces where outside navi- 
gational aids are not available. 

The results presented in Section ITT1 and the above defi- 
nitions lead to the simple conclusion that the Divider set 
is of a curve S in the plane is included in the closure of 
II(S'). Proof of the above statement is provided if a series 
of well known results from differential geometry arc given 
here in the form of lemmas, without proof. 

Lemma III. 3. // the evolute E of a thrice continuously 
differentiable curve S in 5ft 2 is considered, the cusps of the 
evolute correspond to points of minimum or maximum 
curvature of the osculating circles of S at those points. 
The cusps corresponding to maximum curvature have the 
property that if a small open disk B(k, e) is taken with the 
cusp as a center, then the cusp is closest to S than any 
other point of the evolute within the disk. If an analogous 
disk is taken at a cusp corresponding to minimum cur- 
vature, then the cusp is furthest from S than any other 
point of the evolute within the disk. 

Lemma III. 4. Let a curve S in 5ft 2 and its evolute E 
be considered. By the theory developed in Q/ ; n(S') is 
the area in 5ft 2 in the convex side of the evolute. If S is 
homeomorphic to the area ofH(S) will be enclosed 
within the evolute. If, on the other hand, S is a curve 
in 5ft 2 and is homeomorphic to (0, 1), (0, 1] or [0, 1], then 
H(S) is the area of 5ft 2 toward which the evolute E turns 
its convex side. 

A typical example is an ellipse ITTll 

If, on the other hand, S is a curve in 5ft 2 and is home- 
omorphic to (0, 1), (0, 1] or [0, 1], then II(S') is the area 
of 5ft 2 toward which the evolute E turns its convex side. 
A typical example is a parabola IIIII 

Proof. By the theory developed in [| , if S is homeomor- 
phic to a one - dimensional connected subset of a straight 




FIG. 1: The set n(5) of an ellipse S is the set of 
points within the interior of the area within the evo- 
lute, called, in this case an asteroid. It is indicated 
by the red line. It is the set of all points in 5ft 2 for 
which the curvature of locally convex type takes a non 
zero value. The Divider of S is the segment of the 
great axis within the closure of II(S'). Its endpoints are 
the centers of osculating circles of maximum curvature. 
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FIG. 2: The evolute of a parabola is called a Neile's 
parabola. The set H(S) lies within the area delim- 
ited by the red line, where the evolute turns its con- 
vex side. The Divider of a parabola is the part of its 
axis starting from the center of the osculating circle at 
the vertex of maximum curvature and going to infinity. 



line, then a necessary and sufficient condition for the re- 
lation p G II(S') to be valid, is that there is at least one 
line p and q normal to S, such that the distance \p, q\ of 
p from q is a local maximum of the distances of p from 
points of S. If, on the other hand, S is homeomorphic 
to then two line segments \p, qi\, \p, q2 , which must 
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be local maxima of the distance of p from points of 5 |l| . 
For the above to hold, there must exist analogous tan- 
gent lines from p to the evolute E of 5. For the distance 
of p from a point q of 5 to be a local maximum, the point 
of contact t of pq with E must lie between p and q. The 
above statement can be true only if p is at the convex 
side of E, so that the appropriate tangent or tangents 
may exist. □ 

Lemma III. 5. If S is homeomorphic to S^' and dif- 
ferent from S^ 1 ' , then from every point p ofH(S), there 
will be at least four tangents to E from p. Correspond- 
ingly, they will be normal to S, by definition of involutes 
and evolutes. Two of them will define local minimum 
distances of p from 5 and two will define local maximum 
distances from p to S MIR Each one of then is tangent 
to E at a point t{, i = 1, 2, 3, 4, which is the center of an 
osculating circle of S at a point q £ 5 . The lines defin- 
ing local maxima and minima of distance will alternate 
in succession as they meet 5. If, on the other hand, S 
is a curve in 3? 2 and is homeomorphic to (0, 1), (0, 1] or 
[0, 1] but is not a straight line subset, then there will be at 
least three tangents from p to E. Two of them will define 
local minimum distances of p from 5 and one of them 
will define a local maximum distance. The line defining 
a local maximum will lie between the other two lines. JJ/ 

Theorem III. 6. If a tangent is defined from p £ 11(5) to 
the evolute E of S, being normal to 5 at q, then the points 
p' on the straight line (p, q) have the following properties: 
If a point p' lies on the far side of q relative to the cen- 
ter t of the osculating circle, the distance \p, q\ is a local 
minimum, in the sense that if a small enough open neigh- 
borhood of q, N(q), is considered on S, then Vg' £ IT(5) 
the relation: \p,q\ < \p,q'\ holds. 

The same is true if p lies between q and t. If p is at 
the far side of t relative to q, in other words if t is be- 
tween p and q then the distance \p, q\ is a local maximum, 
meaning that, if a small enough open neighborhood of 
q, N(q), is considered on 5, then Mq 1 £ 11(5') the relation: 
\p, q\ > \p, q'\ holds. Furthermore, if p = t, then if t is the 
center of an osculating circle with maximum curvature, 
as described above, the relation \t, q\ < \t, q'\ holds. If t is 
the center of an osculating circle of minimum curvature, 
the relation \t,q\ > \t,q'\ holds. Finally, if t is the center 
of an osculating circle with neither maximum nor min- 
imum curvature, the distance \t, q\ is neither maximum 
nor minimum. If any neighborhood N(q) is considered 
on 5, however small, the points q' in N(q) for which the 
osculating curvature is larger than that in q will lie out- 
side the osculating circle at q, while the points q' where 
the osculating curvature is larger than that of q will lie 
within the osculating circle ay q. The above results are 
well known facts from differential geometry. An osculat- 
ing circle of a curve at a point p is in contact of at least 
the second degree with the curve at p. If the curvature 
is locally maximum or minimum, then the contact is at 
least of the third degree. If the curvature is locally a 



maximum, then the curve lies outside the osculating cir- 
cle at a neighbourhood of p. If the curvature is a local 
minimum, then the curve lies inside the osculating circle 
at a neighbourhood of p. In ordinary cases, i.e., when the 
contact at the point in question does not happen to be 
of an order higher than the second, the circle of curva- 
ture will not merely touch the curve, but will also cross 
it |17| . These results can be naturally extended to curves 
and surfaces in three dimensions 0, Il7| . 
Now an important Theorem will be presented: 

Theorem III. 7. Let S be a curve in K 2 , homeomorphic 
to (0,1), twice continuously differentiable, unbounded, 
not containing arcs of constant curvature and with only 
one point qo of maximum curvature. The Divider of 5 is 
contained in the closure ofH(S), at the part o/5R 2 where 
its evolute E turns its convex side. The end points of the 
Divider are cusps of the evolute E of S , corresponding to 
maximum curvature of the relative osculating circles. 

Proof. If 5 is as above, obviously the evolute E will have 
two branches joining at a cusp t and having a common 
tangent there. The cusp t will be the center of an os- 
culating circle of maximum osculating curvature. If a 
point p ^ t belongs to the Divider, there will be exactly 
three straight line segments \p, q\ | , \p, q 2 1 , \p, qz | , normal 
to 5 and tangent to E. One of them, \p, q 2 \, will define 
a local (and global) maximum distance of p from 5 and 
will lie among the other two. The other two distances, 
\p, qi\, \p,q3\, will be equal by definition of the Divider. 
Also, being equal to the radius of the supremum of disks 
having contact with 5 at only one point, it will also be 
the infimum of the disks having more than one contact 
points with 5. The set B {2 \p, \p, q 2 \) f| 5 will be dis- 
connected, containing only two distinct points qi and q%. 
Therefore, p £ 11(5) On the other hand, let the cusp t of 
E be considered. It will be the center of the circle of max- 
imum osculating curvature of 5 at q. As such, by Lemma 
IIII.4l and IIII.6I it will also be the center of the maximum 
disk contacting 5 at its single point q. Therefore, the 
disk having t as a center and \t, q\ as a radius, will define 
the contact curvature of 5 at q. The cusp t will belong 
to the Divider and since it will not belong to 11(5), it 
will be an end point of the Divider. The osculating disk 
will not intersect 5 at any other point, since 5 has only 
one point of maximum curvature. Similar results can be 
easily proved for more general curves, homeomorphic to 
(0, 1], [0, 1], 5' 1 ', lines which are receptive of a normal at 
each point but their osculating curvature function is piece 
by piece continuous, or even disconnected sets of curves. 
The above results do not hold only in cases where S has 
points of self intersection IIIII □ 



All the above statements can be verified by the Divider 
defining equations presented below: Let a curve 5 in be 
defined by the parametric functions: 
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FIG. 3: A hypotrochoid S and its evolute are the best example of an exception to the rules for curves without double 
points. Note that vertices of a curve correspond to its evolute cusps. The lines are radii of tangent circles. The area 
inside the evolute defines the set n(S') of the hypotrochoid S. The points at the interior have a positive curvature of 
locally convex type. The Divider interior points are contained in the set n(5). The cusps of the evolute corresponding 
to the maximum curvature points of the hypotrochoid, are not end points of the Divider. This is because the oscu- 
lating disks at the maximum curvature points are larger than the corresponding maximal disks in the same points. The 
cross points of the hypotrochoid arc points of the Divider with zero radius maximal disks, having infinite curvature of 
locally convex type. The interior of the smaller curved asteroid decagon inside the hypotrochoid is the area where the 
disks defining the curvature of locally convex type may have an intersection with the hypotrochoid with five distinct 
connected components. As a special case, the central junction point of the Divider, situated at the center of the Figure, 
is the starting point of five branches of the Divider, these being straight lines radiating from the junction to each cross 
point of the hypotrochoid. There are five more branches of the Divider, uniting the hypotrochoid cross points within 
each outer lobe of the hypotrochoid. The last five branches of the Divider are straight lines starting from each cross 
point of the hypotrochoid, passing through the corresponding cusp of the evolute which is a center of a minimum curva- 
ture osculating disk and going on to infinity. In this case the Divider is a multiply connected graph containing singular 
points of zero contact curvature due to the existence of analogous singularities (cross points) in the hypotrochoid. 



x\ 

X 2 



x 2 {t) 



(5) 
(G) 



Given any value t\ of the defining parameter of 5, 
the following quantities will be calculated, by the equa- 
tions and inequalities presented below: x\{t\), x 2 {ti), t 2 , 

X\(t2), X 2 (t 2 ), Xio, £20- 



x\ = xi(h) 

xi = x 2 (h) 

x\ = Xi(t 2 ) 

x 2 = x 2 {t 2 ) 



(7) 
(8) 
(9) 
(10) 



/ , ,dxi(ti) , , , .dx 2 (ti) 
{x 1 {t 1 )-x 10 )—^ + {x 2 (t 1 )-X2o)—£ 11 = 
ati ati 



(xi(t 2 ) - Xio) 



dxi(t 2 
dt 2 



(x 2 (t 2 



x 2 o) 



^dx 2 (t 2 
dt 2 



(11) 
(12) 
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(xi(ti) - x w ) 2 + (x 2 (*i) - x 20 ) = (xi(t 2 ) - x w ) 2 + (x 2 (t 2 ) - x 2 oY 



(13) 



< (xi(ti) - x w ) 



< (x 1 (t 2 ) - x w ) 



d 2 x x {h) , f d Xl (h) 



dt\ 



dt\ 



d 2 xi(t 2 ) ( dxi(t 2 ) 



+ {x 2 (h) - x 20 ) 



+ (x 2 (t 2 ) - x 20 ) 



d 2 x 2 {ti) i d x 2 (tj) 
dh 



dt\ 



dt 2 



d 2 x 2 (t 2 ) ( dx 2 (t 2 ) 



V *2 



(14) 
(15) 



There are important cases where the equations: 



(xi(<i) - x w ) 



(xi(t 2 ) - x w ) 



d 2 x 1 (ti) ( dx\(ti) 



dtl 
d 2 xi(t 2 



dh 

dxx(t 2 ) 
dh 



+ (x 2 {h) - x 20 ) 



+ (x 2 (t 2 ) - x 20 ) 



d 2 x 2 (h 
dtj 

d 2 x 2 (t 2 
dt 2 . 



dx 2 (ti 
dh 

dx 2 (t 2 
dh 



= 



(16) 
(17) 



r 



are true and yet the distances \p, q\ | , \p, q 2 | are local 
minima. In that case, the sign and order of the lowest, in 
order succession, nonzero derivative decides if \p, q\ | and 
\p, q 2 1 are local maxima or local minima of the distances 
of p from the points q of 5*0. This simplified version of 
minimum distance conditions is given here for reasons of 
space. In cases of curves or surfaces with boundaries, the 
existence of local maxima or minima of the distance of p 
from S is not necessarily connected with the existence of 
normals from p to S. In such cases, minimum or maxi- 
mum distances may lie on the boundary of S, in which 
case a set of relations analogous to ^ [21 131 Q may not 
exist, or may yield only partial results, in the creation 
of the Divider. In such cases, the algorithms for the cre- 
ation of the Divider should use procedures not entirely 
based on such well defined and solvable equations. 

As h spans S, the coordinates xio, x 2 q, of the Divider 
are calculated, as functions of the variable t\. 

Let a curve S in 5ft 2 , homeomorphic but not similar to 
S{\), be considered. The above equations are consistent 
with the definition of maximal disks 0, 0, B IS EH 113 • 
The maximal disks definition and the Divider definition 
are equivalent as long as the Euclidean metric ^ is used. 
The definition of contact disks and maximal disks as 
stated in this work and in are not equivalent to the 
traditional definition of maximal disks, if one of the other 
metrics, |21 and |3| is used. 

The above equations 171 to IT31 arc designed to find the 
coordinates xio, of all points x\(b\) and x 2 {h) which be- 
long to the Divider of S. The procedure is as follows. If a 
value h of the parameter t in the relations \J\ and [H] defin- 
ing S is given, then, by [7] and 00 a point q\ with coordi- 
nates x\(h) and x 2 {h) will be defined. Equations and 
1101 define a new value t 2 of t and the coordinates xi(t 2 ) 
and £2(^2) °f a corresponding point q 2l with appropriate 
properties designated by the following equations. The 



rest of the equations, ^| to ^3 define the coordinates 
xqi and Xq 2 of a point p £ S 2 which, as stated above, 
will belong to the Divider of S. Eg uat ion 1111 signifies that 
\p, qi \ is normal to S. Eauation ll2l signifies that \p, q 2 \ is 
normal to S. Equation 1131 signifies that \p,qi\ = \p, q 2 \- 
Finally, inequalities ^] ^| certify that \p, qi \ and \p, q 2 \ 
are local minima of the distances of p from the points of 
S. 

The above equations can be easily extended to curves 
and surfaces in 5ft 3 . Other, more general geometric ob- 
jects in abstract metric spaces may be defined by similar 
equations. By the theory of curves and surfaces, if S, 
be it a curve or a surface in 5ft 2 or 5ft 3 , is unbounded in 
the topological sense and, furthermore, has no boundary 
points by its intrinsic topology [Tsj . the set of equations 
[7J to E| has at least one solution. In 5ft 2 , this would mean 
a minimum of distances \p,q\, q G S, if one of the in- 
equalities El or El would hold. Conversely, if none of 
the above inequalities holds and the solution yields a lo- 
cal maximum, there will be at least two more solutions, 
both defining local minima of distances \p,q\, q G S, at 
points qi, q 2 . By El \PiQi\ = [Pjffel- Then the point p 
would belong to the Divider of S. Analogous results hold 
for bounded curves or surfaces, homeomorphic to S(l) or 
5(2) but not having everywhere constant curvature. A 
special case is that of a center c of an osculating circle 
of constant curvature at every point q' on N(q). In that 
case the following are true: 

1. S contains a compact, connected circular arc Q with 
c as its center. 

2. Every contact direction for every q G Q passes 
through c. 

3. The disk of the osculating circle is the supremum of 
the disks which have with S only one common point 
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q, belonging to Q. The center c of the osculating 
circle for the points q of Q is the corresponding 
point for the Divider for Q. 

All the above results are presented here without proof, 
since they are well known from classical differential 
geometry 0- 

IV. DISCUSSION AND CONCLUSIONS 

In the above, the main contribution is the definition 
of a new skeleton concept, the Divider set. Being in 
some sense the reverse of the maximal disk definition, 
the definition of the Divider is given in reference to the 
points of a curve or surface. This may or may not be 
the boundary of an area, finite or infinite. The defini- 
tion and the construction of a maximal contact disk will 
start from a point on the curve or surface. On the other 
hand, the classical definition of Blum, has to do with 
maximal inscribed disks. Therefore, Therefore the algo- 
rithms creating it will be in principle more efficient from 
algorithms referring to all points inside a given enclosed 
area, as are the grassfire or wave front algorithms and 
most maximal disks algorithms these authors are aware 
of@0,ESIIl[n],[llEm- This seems to be the 
case in the first attempts for a Divider creating algorithm 
by the use of the defining equations using the Euclidean 
metric 11111 Mil Mil and also in the case of the discrete 
lattice |5j. In that last case, the maximum coordinate 
metric is used. In most cases in the literature, for a small 
sample see [H [H [3 , the Euclidean metric is used for 
the creation of skeletons, even in a discrete lattice envi- 
ronment. 

The maximal disk definition, as mentioned above, is 
not well suited for the maximum coordinate metric. On 
the other hand the Divider definition can be naturally 
modified [f| for a discrete lattice environment. The re- 
sult is a connected set of cells, having at most a two cells 



width wherever the distance of the lines is an even num- 
ber of cells. This is a well known problem with discrete 
lattices 0, 13 and there are many issues to be dis- 
cussed. One of them, arguably the most important, is 
the creation of a skeleton suitable for image compression 
and more or less faithful restoration 0,0, IS 01- The au- 
thors of this work are of the opinion that in some cases 
reproduction of a faithful image requires that all skeleton 
pixels should be kept and no further thinning algorithm 
should be applied. In other cases, like OCR, where no 
faithful reconstruction of the image is required but the 
objective is the extraction and preservation of some im- 
portant features of the initial image, there are some sim- 
ple thinning algorithms that can be applied with good 
results 0- 

In general, the Divider seems to work in a satisfac- 
tory manner for the applications of interest to the au- 
thors, for example OCR and robotic navigation in en- 
closed environments, among others. To conclude, by the 
work so far, the Divider concept seems to have some 
distinct advantages compared with many other state of 
the art skeletonization methods. It has a precise mathe- 
matical definition, easily implemented algorithms for two 
and three dimensional Euclidean spaces utilizing the Eu- 
clidean metric, as well as discrete lattice spaces, utilizing 
the maximum coordinate metric. It seems to be promis- 
ing for many applications, although further testing and 
comparison with alternate methods is still to be done 
in future work. It has especially promising attributes in 
specific interesting applications, such as handwritten text 
recognition and robotic navigation, among others. 
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[1] If the set of such radii r for which the above relation holds 
is empty, which means that there are no closed balls such 
that their intersection with S will be disconnected, then 

K Lc.t = °. as the 

inverse of the infimum of an empty set. 

The simplest examples of such a situation are plane or 

spherical surfaces. 
[2] In the definition is given for manifolds M d in W 1 . It 

can be applied equally well to curves S in 5R 2 or K 3 , as 

well as surfaces in 5ft 3 . 
[3] For a detailed proof, see [J. 
[4] Proof of all above results in 0. 

[5] By definition, every tangent to the evolute E will be nor- 
mal to S at q and the osculating curvature of S at q will 
be k(q) = l/|t, g|. A theorem, including results well known 
from the theory of differential geometry of curves and 
surfaces, will be presented below. Similar results are pre- 
sented and proved in the relative literature (For example, 
see Ili). 

[6] Relative theorems are contained in most reference books 

of infinitesimal calculus 
[7] See also pertinent results in 0,0 



